Low density spin-polarized transport in 2D semiconductor structures: The enigma of 
temperature dependent magnetoresistance of Si MOSFETs in an in-plane applied 

magnetic field 

S. Das Sarma and E. H. Hwang 
Condensed Matter Theory Center, Department of Physics, University of Maryland, College Park, Maryland 20742-4111 

(February 2, 2008) 

The temperature dependence of 2D magnetoresistance in an applied in-plane magnetic field is theo- 
retically considered for electrons in Si MOSFETs within the screening theory for long-range charged 
impurity scattering limited carrier transport. In agreement with recent experimental observations 
we find an essentially temperature independent magnetoresistivity for carrier densities well into the 
2D metallic regime due to the field-induced lifting of spin and, perhaps, valley degeneracies. In 
particular the metallic temperature dependence of the ballistic magnetoresistance is strongly sup- 
pressed around the zero-temperature critical magnetic field (Ba) for full spin-polarization, with the 
metallic temperature dependence strongest at B = 0, weakest around B ^ Bs, and intermediate at 
B > Bs. 
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I. INTRODUCTION 

The phenomena of 2D metallic behavior and the as- 
sociated 2D metal-insulator transition (MIT) comprise 
a complex set of low-temperature transport behavior 
of low-density 2D electron (or hole) systems in high- 
quality (i.e. low-disorder and high-mobility) semicon- 
ductor structures. In particular, the 2D resistivity, 
p{n,T, B), shows intriguing and anomalous dependence 
on carrier density (n), temperature (T), and applied mag- 
netic field {B) parallel to the 2D plane [1,2] . For example, 
the low temperature p{T), in zero applied field, shows 
remarkably strong "metallic" -like (i.e. dp/dT > for 
n > Tie) temperature dependence for 2D carrier densi- 
ties above the so-called critical density (ric) for the 2D 
MIT whereas, for n < Uc, the system exhibits insulating 
behavior {dp/dT < 0). The strength of the metallicity 
(i.e. how strongly p{T) increases with T at low temper- 
atures for n > He) and the actual value of ric are highly 
nonuniversal, and depend strongly (and nontrivially) on 
the material and on the 2D sample quality. The applica- 
tion of an in-planc magnetic field B has several interest- 
ing effects on the 2D metallic phase [3-10]: (1) at a fixed 
low T, the system develops a large positive magnetore- 
sistance with p{B) increasing very strongly (by as much 
as a factor of 4) with B upto a maximum field Bg] (2) for 
B > Bs, p{B) either saturates (or increases slowly with 
B for B > Bs) showing a distinct kink at B = Bs; (3) 
the temperature dependence of resistivity is considerably 
suppressed in a finite parallel field; (4) for 2D electrons 
in Si MOSFETs, which is the most extensively studied 
system in the context of 2D MIT, the external magnetic 
field at densities close to the critical density (i.e. n > ric) 
seems to drive the system into a strongly insulating phase 
generating a huge positive magnetoresistance and mani- 
festing the strongly insulating activated temperature de- 
pendence typical of the zero- field transport for n < ric- 



In this paper we restrict ourselves to the highly metal- 
lic ( "ballistic" ) n > 2nc regime and investigate theoreti- 
cally the enigmatic temperature dependence of the paral- 
lel field magnetoresistatce p{B; T) of n-Si MOSFETs as 
described in items (1) ~(3) above. The item (4), which 
arises from the destruction of the 2D effective metallic 
phase by a large in-plane field in the n > ric regime, is 
obviously related to the 2D metal-insulator transition it- 
self (or more precisely, its parallel field dependence) and 
is beyond the scope of the current theory which deals 
only with the apparent effective 2D metallic phase in the 
ballistic regime. We include in the theory only resistive 
carrier scattering by screened effective disorder arising 
from charged impurity centers randomly distributed at 
the Si-Si02 interface with the screening calculated within 
the finite wave vector self-consistent field random-phase- 
approximation (RPA) at finite temperatures fully incor- 
porating the carrier spin polarization induced by the ap- 
plied in-plane magnetic field. In our theory the magnetic 
field effects enter only through the carrier spin polariza- 
tion correction. The temperature dependence (at zero 
field) [11] and the parallel-field dependence (at T = 0) 
[12] of p{T; B) have earlier been individually theoreti- 
cally studied in the literature for Si MOSFETs within 
the screening theory formalism - our goal here is to 
combine the two to obtain a full theory for p{T] B) in 
Si MOSFETs including spin-polarization effects on the 
finite-temperature screening properties. 

We emphasize that, although the 2D "metallic" trans- 
port properties are anomalous in the presence of a par- 
allel magnetic field in all of the 2D systems studied so 
far, the 2D Si MOS system is unique in exhibiting a par- 
ticularly intriguing temperature dependent magnetore- 
sistance p{T,B) where p{T) for B > Bs and n > l.bric 
seems to be essentially a constant at low temperatures 
[7,8]. In this paper we propose a simple physically moti- 
vated theoretical explanation for the anomalous trans- 
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port properties of Si MOSFETs in the presence of a 
strong parallel magnetic field. Our explanation is based 
on the substantial modification of the effective quenched 
disorder in the system (arising from the screened back- 
ground impurity potential) due to the applied paral- 
lel field which strongly affects the 2D screening prop- 
erties. We believe (and show in this paper) that a qual- 
itative understanding of the temperature dependence of 
p{T; B) can be obtained on the basis of an effective field- 
dependent disorder. 

The applied parallel magnetic field has, in principle, 
seven distinct (and sometimes, opposing) physical effects 
on the 2D metallic transport behavior: (1) the magneto- 
orbital coupling [13] due to the finite thickness of the 2D 
layer leads to an anisotropically increasing 2D effective 
mass, consequently giving rise to a positive anisotropic 
magnetoresistance; (2) the parallel field couples the 2D 
in-plane dynamics to the dynamics in the confinement 
direction (perpendicular to the 2D plane) leading to in- 
tersubband scattering (i.e. 2D to 3D crossover) induced 
positive magnetoresistance [13]; (3) the parallel field may 
effectively "enhance" weak localization "insulating" tem- 
perature correction [14] by suppressing various "metal- 
lic" contributions (e.g. screening) to the temperature de- 
pendent resistivity; (4) the parallel field-induced carrier 
spin polarization leads to an enhancement of the effec- 
tive Fermi momentum, kpiB), which tends to contribute 
a negative magnetoresistance through the decrease of the 
charged impurity scattering matrix elements; (5) the par- 
allel field-induced carrier spin polarization reduces the 
strength of screening (by as much as a factor of 2) since 
the electronic density of states decreases from 2 to 1 as 
the 2D electrons become completely spin-polarized; (6) 
the parallel field may further suppress screening (by as 
much as an additional factor of 2) in Si MOSFETs by 
lifting the valley degeneracy factor (for example, from 
2 to 1 if Ay > Ef, ksT, where A^, is the valley split- 
ting in the presence of the parallel field); (7) the parallel 
field induced modification of the electron-electron inter- 
action [15], due to the spin-polarization of the 2D elec- 
trons which suppresses the "triplet" interaction channel. 

We develop a theory for piT, B, n) in Si MOSFETs in- 
cluding the three physical mechanisms [(4) - (6) listed 
above] which are important for Si MOS systems. We 
leave out the magneto-orbital effects (items (1) and (2) 
above) [13] because the magneto-orbital corrections are 
rather small in Si MOSFETs since the quasi-2D layer 
width (i.e. the confinement size transverse to the 2D 
plane) is rather small, making magneto-orbital coupling 
essentially negligible in Si MOSFETs except at very high 
parallel fields. We do not consider the weak localization 
correction since it is straightforward to include it (in an 
ad hoc fashion) if experimental results warrant such a 
theoretical adjustment [14]. We also uncritically ignore 
all electron-electron interaction corrections beyond the 
screening (i.e. ring diagrams) effects, arguing screening 



to be the dominant physical mechanism controlling trans- 
port limited by Coulomb disorder arising from charged 
impurity scattering. Interaction effects in the presence 
of a parallel magnetic field have been considered in the 
literature [15]. 

The main physical effect we consider, namely the par- 
allel field induced decrease of the 2D density of states 
leading to a strong field induced suppression of screening 
(and hence a strong positive magnetoresistance), has ear- 
lier been discussed in the literature [12], primarily in the 
context of the T = magnetoresistance itself (whereas 
our focus is on the temperature dependence of magne- 
toresistance). We have recently [16] discussed the similar- 
ity between the behavior of p{B; T = 0) and p{T; B = 0) 
in Si MOSFETs [17] as arising from the field-induced or 
the temperature-induced suppression of screening, lead- 
ing to qualitatively "similar" parallel-field dependent (at 
T = 0) and temperature-dependent (at B = 0) effective 
disorder in MOSFETs. (As an aside we mention that the 
situation is quite different in 2D GaAs electron [18] and 
hole [19] systems where substantial magneto-orbital cou- 
pling is present, leading to qualitatively different p{B) 
and p{T) behaviors.) The goal of the current paper is 
to develop a theory for p{T;B), the temperature depen- 
dence of the 2D magnetoresistivity in Si MOSFETs, par- 
ticularly at large fields B > Bg, where the 2D system is 
presumably completely spin-polarized. The motivation 
for our theoretical work comes from the recent experi- 
mental work [7] reporting a puzzling absence of any tem- 
perature dependence in p{T;B > B^) of Si MOSFETs 
in the completely spin-polarized large applied parallel 
field {B > Bg) ballistic (n > 1.5nc) transport regime. 
Since the zero-field (i.e., spin unpolarized) 2D transport 
in these high-quality low-density Si MOSFETs is charac- 
terized by strong metallicity, i.e. a strong metallic tem- 
perature dependence in p(T), the apparent complete sup- 
pression of the metallicity in the spin-polarized system is 
intriguing and has attracted a great deal of attention. 
The complete reported absence of any temperature de- 
pendence of p(r; B > Bs) in Si MOSFETs [7] also stands 
in sharp contrast to the corresponding situation [20] in n- 
Si/SiGe 2D structures, where p(T; B > Bg) shows metal- 
lic temperature dependence in the spin-polarized case, 
albeit with a reduced magnitude of dp/dT in accordance 
with the screening theory. We discuss these issues in 
great details in the next four sections of this article. In 
sec. II we provide the Boltzmann transport theory. In 
sec. Ill we give our calculated results. In sec. IV we 
discuss our results. We conclude our paper in sec. V. 

II. THEORY 

We use the semiclassical Boltzmann transport theory 
[11] including only the effect of resistive scattering by 
random charged impurities at the Si-Si02 interface - 
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at low carrier densities (and at low temperatures where 
phonon scattering is unimportant) of interest in the 2D 
MIT problem, transport in high-mobility Si MOSFETs 
is known to be predominantly (but perhaps not entirely) 
limited by long-range charged impurity scattering. The 
density of the random interface charged impurity centers 
is therefore the only unknown parameter in our model, 
which sets the overall scale of p{T = 0) without affect- 
ing the (T, n, B) dependence of resistivity that is of in- 
terest in the problem. The bare disorder potential aris- 
ing from oxide charged impurity centers being long-range 
Coulombic in nature, the most important physics ingredi- 
ent (at the zeroth-order) that any transport theory must 
incorporate is the regularization of the bare long-range 
Coulombic disorder potential by screening the impurity 
potential. Within our zeroth-order minimal transport 
theory, this is precisely what we do by calculating the fi- 
nite temperature screened effective disorder through the 
static (finite-temperature and finite-wave vector) RPA 
screening of the bare long-range Coulombic disorder [11]. 
We then use the Boltzmann theory to calculate the finite- 
temperature carrier resistivity limited by scattering due 
to the effective "regularized" (i.e. screened) impurity dis- 
order. We include in our theory the realistic effects of the 
quasi-2D layer width of Si MOSFETs through the appro- 
priate quantum form factors, and most importantly we 
include the nonperturbative effect of the parallel mag- 
netic field through the spin polarization of the 2D elec- 
tron system fully incorporating the physical mechanisms 
(4) and (5) discussed above [12,16] in Sec. I. (The in- 
clusion of the mechanism (6) in the theory, namely the 
field-induced lifting of valley degeneracy, i.e. valley polar- 
ization, is straightforward by adjusting the valley degen- 
eracy factor gv in the 2D density of states and is discussed 
below.) 

The formal aspects of the calculation for p{n; T; B) in 
the presence of the parallel magnetic field are essentially 
the same as those for p{n] T) at B = except for the 
existence of different spin-polarized carrier densities n±. 
When the parallel magnetic field is applied to the system 
the carrier densities n± for spin up/down are not equal. 
Note that the total density n = n+ + ri- is fixed. The 
spin-polarized densities themselves are obtained from the 
relative shifts (i.e. the spin-splitting) in the spin up and 
down bands introduced by the Zeeman splitting asso- 
ciated with the external applied field B. At T = 0, 
one simply has n± = n(l ± B/Bs)/2 for B < Bg with 
n+ = n and n_ = for B > Bg (and n+ = = n/2 at 
B = 0), where B^, the so-called saturation (or the spin- 
polarization) field, is given by gfisBg = 2Ep where g 
is the electron spin g-factor and pB the Bohr magneton. 
For T ^ 0, n± is determined using the Fermi distribution 
function in the standard manner. Thus, in the presence 
of the magnetic field the total conductivity can be ex- 



pressed as a sum of spin up/down contributions 

a = a++a-, (1) 

where a± is the conductivity of the (±) spin subband. 

The total carrier resistivity p is defined by p = 1 /a. The 
conductivities a± are given by 



a± = 



n±e'^{T±) 



(2) 



where m is the carrier effective mass, and the energy aver- 
aged transport relaxation time (t±) for the (±) subbands 
are given in the Boltzmann theory by 



{r±) = 



J d£T{e)e 
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de 


J dee 
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(3) 



where r(e) is the energy dependent relaxation time, and 
/±(e) is the carrier (Fermi) distribution function 



(4) 



where /3 = {kBT)~^, and p± = ^ \n[-l + exp{i3Ep±)] 
is the chemical potential for the up/down spin state 
(with Fermi energy Ep± ) at finite temperature. The 
spin-polarized transport can then be calculated within 
this two-component (spin-up and -down) carrier system 
(without any spin-fiip scattering since the bare impurities 
are non-magnetic), with the screening of the bare disor- 
der being provided by both spin up and down carriers. 

Within the framework of linear transport theory the 
relaxation time for a carrier with energy £fe is given by 



1 27r 



/OC' 
dzN^''\z)\u^"\k--k'-z)\'' 
a,K- -°° 

X(l - C0s6'kk')(5(£fc - Efc/), (5) 



where ^kk' is the scattering angle between wave vectors k 
and k' and Sk = h^k'^/2m; the screened scattering poten- 
tial is denoted by u^"\q; z) with q = k — k' and z is the 
confinement direction normal to the 2D layer. N^^\z) 
in Eq. (5) is the 3D charged impurity density of the a-th 
kind of charged center. Here we have assumed that the 
charged centers are distributed completely at random in 
the Si-Si02 interface for MOSFETs. The screened impu- 
rity potential u'^"^ (q; z) is given by: 



27rZWe^ 



Ft>{q;z), 



(6) 



where Z'"^ is the impurity charge strength, R is the 
background (static) lattice dielectric constant, and Fi is 
a form factor determined by the location of the impu- 
rity and the subband wavefunction tl}{z) defining the 2D 
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confinement. The finite wave vector dielectric screening 
function e{q) is written in the RPA as 



27re2 

e(g) = l f{qMq,T), 



(7) 



where f{q) is the Coulomb form factor arising from the 
subband wavefunctions tp{z). In the strict 2D limit 
f{q) = 1. Il{q, T) is the 2D irreducible finite-temperature 
(and finite wave vector) polarizability function and is give 
by n(9;T) = n+(g;T) +n_(g;T), where n±(g;T) are 
the polarizabilitics of the polarized up/down spin states 
(±). At finite temperature wc have 



n°±(g,A^') 



cosh^ |(/z± - fi') ' 



(8) 



where Ilj.{q,Ep) = W±{q) is the zero-temperature non- 
interacting static polarizability, given by 



Ul{q) = NF 



1 - ^Jl-{2kF±/qfO{q - 2kF±] 



(9) 



where Np = g^m/2Tr is the density of states per spin at 
Fermi energy, and kp± = (47rn±/gf„)^/^ is the 2D Fermi 
wave vector for the spin up/down carriers. Note that 
Qy (=2 in the usual B = Si MOS case) is the val- 
ley degeneracy, and the spin degeneracy, by definition, is 
assumed to be lifted by the in-plane field B, the usual 
unpolarized B = paramagnetic case being kp+ = kp- ; 
n+ = n_ = n/2. 

We mention that Eqs. (1) (9) for the finite- 
temperature carrier transport properties in 2D systems 
comprise a complex set of multi-dimensional integrals 
along with the determination of the quasi-2D Coulomb 
form factors as well as the chemical potential of the 
system. In the asymptotic regimes of T/Tp <^ 1 and 
T/Tp > 1, Eqs. (1) (7) simplify (as discussed in 
rcf. [11]) giving rise to simple analytic behaviors in 
piT) ~ 0{T/Tp) for T/Tp «; 1 and p{T) ~ 0{Tp/T) for 
T/Tp ^ 1, but these asymptotic analytic behaviors are 
of limited experimental relevance since few experimental 
parameter regimes satisfy the required conditions [11] for 
the asymptotic behavior. Also, at the lowest experimen- 
tal temperatures p{T) invariably saturates essentially in 
all experiments. For arbitrary T and n (as well as B) one 
must evaluate Eqs. (1) (7) with sufficient accuracy to 
obtain p{T^ n, B) for 2D systems. 

The most salient aspects of the parallel field induced 
carrier spin polarization are an enhancement of the 2D 
Fermi wave vector by a factor of ^/2 and a suppression 
of the 2D screening ( "Thomas- Fermi" ) wave vector by 
a factor of 2 as the unpolarized ("paramagnetic") 2D 
system becomes completely spin-polarized ("ferromag- 
netic") with the parallel field increasing from i? = to 
B > Bg (with Bs being the full or saturation spin polar- 
ization field) . In the most naive theoretical level one can 



write (at T = 0) p oc {qrp + 2kp)~'^ in the strong screen- 
ing limit for scattering by screened charged impurities, 
leading to 



p{B > B,) 
P(0) 



< 



gTF(O) + 2kp{0) 

qTF{Bs) + 2kF{B,) 



(10) 



In the usual range of experimental 2D densities qxp ^ 
2kp ("strong screening") in Si MOSFETs, and there- 
fore p{B > B,)/p{0) < [qTFiO)/qTF{Bs)]^ < 4. (As an 
aside, we note that in the opposite limit of weak screen- 
ing, qTP <C 2kp, which may be approximately the sit- 
uation in the high-density 2D n-GaAs system, one gets 
p{Bs)/p{0) < 1/2.) At finite temperatures and with finite 
wave vector screening p{B > Bg)/p{0) is expected to be 
less than 4 in Si MOSFETs as observed experimentally. 
Note that the 2D strong screening condition, qrp ^ 2kp, 
occurs at low carrier densities since kp oc y/n and qxp is 
independent of density in the lowest order. 

The suppression of screening due to the parallel field- 
induced spin polarization has direct consequences for 
the temperature dependence of the magnetoresistivity 
p(r, B) since the temperature dependence of screening 
gives rise to an effective metallic behavior of p{T) at 
zero magnetic field [11,16]. In particular, the strength 
of the metallicity, i.e. how strong the low temperature 
metallic temperature dependence is in a 2D system, de- 
pends on the dimensionless parameter {qTp/2kpY ~ 
{A.2/n)[gsgi,Y in Si MOSFETs where n is the carrier 
density n measured in the units of 10^^ cm~^ and g^, g^, 
are respectively the spin and valley degeneracy factors 
(with gs = 2 and g^ = 2 being the usual zero-magnetic 
field spin and valley mipolarized case). As the paral- 
lel field increases, < B < Bg, the 2D carriers become 
spin-polarized with gg decreasing from 2 (for B = 0) 
to 1 {B > Bg), consequently suppressing the metallic- 
ity parameter [qxp /2kpY by a factor 8 between 5 = 
and B > Bg. This implies that the temperature depen- 
dence of the finite field resistivity p{T,B,n) for B > Bg 
at a particular carrier density Ub will be approximately 
(and qualitatively) similar to the corresponding zero-field 
metallic temperature dependence at a carrier density hq 
which is roughly 8 times higher: no « Sub- Since the 2D 
metallic temperature dependence is strongly suppressed 
by increasing density, these simple screening consider- 
ations immediately suggest a very strong suppression 
of 2D "metallicity" (i.e. the temperature dependence 
of 2D resistivity) at high parallel fields. An equivalent 
physical way of describing this strong screening induced 
(through the spin-polarization dependence of screening) 
suppression of the temperature dependence of 2D parallel 
field magnetoresistivity is to observe that the metallic- 
ity parameter in the 2D n-GaAs system (where = 1) 
is much smaller, [qpp /2kp)'^ ^ {0.9/h)gg, due to the 
much smaller electronic carrier effective mass in GaAs 
{mGaAs = 0.067me and msi = 0.19me), and therefore in 
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FIG. 1. The calculated resistivity of Si MOSFET systems 
as a function of the temperature for various densities, n = 1.0, 
1.25, 1.5, 1.75, 2.0xl0"cm"^ (top to bottom), (a) for g., = 2, 

= 2, (b) for gs = 1, g^ = 2, and (c) for g^ = 1, g„ = 1. 



the presence of a strong parallel field {B > Bg so that 
Qs = 1) the Si MOSFET system at a particular carrier 
density nsi has a "mctallicity" which is roughly equal 
to the corresponding zero-field (i.e. gs = 2) metallicity 
in the 2D n-GaAS system at a density ncaAs ~ 2nsi- 
Since the observed metallic behavior (i.e. the temper- 
ature dependence of p) in 2D n-GaAs system is ex- 
tremely weak except [21] at very low electron densities 
(below 10^°cm^^), the temperature dependence of the 
magnetoresistivity f)(T, B) in Si MOSFETs would thus 
be strongly suppressed at large parallel fields. Wc em- 
phasize that the same physical mechanism, namely the 
suppression of screening due to carrier spin-polarization, 
leading to the strong positive magnotorcsistancc p{B) at 
a fixed low temperature also leads to the strong suppres- 
sion of the temperature dependence of p{B, T) at a finite 
parallel field. 



III. RESULTS 

In Fig. 1 we show our calculated temperature depen- 
dence of n-Si MOS 2D resistivity, p{T), for several car- 
rier densities in the unpolarized zero-field case {gs = 2, 
g^ = 2), Fig. 1(a) as weh as the high field {B > Bs) 
fully spin polarized case (g, = 1, 5i/ = 2), Fig. 1(b). 
(All magneto-orbital effects [13] have been ignored in 
these calculations.) A comparison of Figs. 1(a) and (b) 
immediately demonstrates the strong suppression of 2D 
metallicity (i.e. the temperature dependence of p) in the 
high-field spin-polarized system, particularly at higher 
(n > lO^-'^CTO"^) carrier densities. At lower 2D densi- 
ties, however, our calculated p{T) in the spin-polarized 
system seems to manifest stronger metallicity than that 
observed recently by Tsui et al. [7] and by Shashkin et 
al. [8] who have measured the temperature dependence 
of 2D magnetoresistivity in low-density Si MOSFETs, 
reporting essentially no temperature dependence in the 
spin polarized Si MOS systems. 

One possible reason for this weaker experimental tem- 
perature dependence of p{T, B) could be that the strong 
parallel field lifts the valley dependency (i.e., = 1 in 
the high field situation) as well as the spin degeneracy. 
The long-standing valley degeneracy problem in Si MOS- 
FETs is poorly understood theoretically except that it is 
experimentally well-established that the valley degener- 
acy is typically lifted by a small valley splitting, similar to 
spin splitting, of poorly understood theoretical origin. 
It is also established experimentally that A„ increases 
with decreasing carrier density and increasing magnetic 
field [22]. We speculate that it is possible that at low 
carrier densities in the presence of the parallel magnetic 
field the valley degeneracy is lifted (A^, > Ep) so that 
the Si MOS system becomes a spin and valley polarized 
system (gr^ = = 1) for B » Bg. This is not unreason- 
able since the suppression of screening at large parallel 
fields would lead to strong many-body enhancement of 
valley splitting leading perhaps to the lifting of valley 
degeneracy at low densities and high fields. (Whether 
the valley degeneracy is indeed lifted in the experimen- 
tal MOSFETs in the presence of a strong in-plane field 
B ^ Bs can only be ascertained experimentally - we are 
here suggesting only the theoretical possibility.) In Fig. 
1(c) we show our calculated high field Si MOS piT) in the 
spin and valley polarized situation, to be compared with 
the spin-polarized (but valley-unpolarized) case in Fig. 
1(b) and both spin and valley unpolarized (i.e. B = Q) 
case in Fig. 1(a). The spin and valley polarized theo- 
retical results in Fig. 1(c) are remarkably similar to the 
recent experimental results of Tsui et al. [7] - in fact, our 
results in Fig. 1(c) even reproduce the experimental ob- 
servation, as noted in ref. [7], of a weak negative (i.e. p{T) 
decreasing with increasing T) temperature dependence of 
p{T) arising from the rather small value of qTE/'^^kp is 
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FIG. 2. Calculated p{Bs) / p{0) as a function of density (a) 
for quasi-2D system and (b) for pure 2D system. Here, at 
B = 0, Qs = Qv = 2. Note at B = Ss the spins axe completely 
polarized. The lines corresponding to = 1 (2) indicate that 
the valley degeneracy is (not) lifted at B = Bg. 

the full spin and valley polarized situation which leads to 

a sub- leading negative temperature contribution to p{T). 

In Fig. 2 we compare p(T) at different spin and 
valley degeneracies {gs = = 2; = 1, g^, = 2; 
gs = gv = 1) at different densities, with the observa- 
tion that our <?« = 5^ = 1 results are in good agreement 
with the recent experimental Si MOS data in the high 
field {B > Bg) situation [7,8]. Our speculation of the Si 
valley degeneracy being lifted ((/^ = 1), in addition to the 
spin degeneracy, for B > Bg could be further tested by 
considering p{B) at a fixed temperature, as shown in Fig. 
2. At T = and in the unrealistic and incorrect strictly 
2D limit the lifting of both spin and valley degenera- 
cies will lead to p{Bs)/p{0) ^ 16 as n ^ [Fig. 2(b)]. 
This will be in direct disagreement with experimental 
observations where p{Bs)/p{0) < 4 in Si MOS systems. 
But, in the realistic quasi-2D systems (and at finite tem- 
peratures) we find [Fig. 2(a)] that p{Bg)/ p{Q) < 4 for 
n > 2.0 X 10"cm-2 - for n < 2.0 x lO^^cm-^, we find 
p{Bs)/p{0) > 4 in Si MOS systems in the valley polarized 
situation. 

There is, in fact, a qualitative physical explanation [23] 
for the observed temperature independence of the high- 
field magnetoresistance in Si MOSFETs, which is generic 
and universal in nature and does not invoke the ad hoc 
explanation of the field -induced lifting of the silicon val- 




FIG. 3. Spin-polarization q = (n+ — n-)/(n+ + n-) as a 
function of in-plano magnetic field for n = 1.5 x lO^^crn^^ 
and for different temperatures T = 0, 1.0, 2.0, 3. OK. Bs is the 
T=OK full spin-polarization field. Note that for T 0, full 
spin-polarization requires B w 2Bs. 

ley degeneracy we propose above as a possibility. This 
explanation is, however, based on a categorical repudia- 
tion of the existing experimental claims [7] that the con- 
stancy (i.e. temperature independence) of p{T; B > Bg) 
as a function of temperature at high parallel fields does 
not in any way indicate a fundamental suppression of 
"metallicity" (i.e. the metallic temperature dependence) 
in the fully spin-polarized Si MOSFETs, as has been re- 
peatedly emphasized by several experimental groups in 
the past [7,8,1]. Instead, the rather generic explanation, 
recently suggested in ref. [23], is that the observed (essen- 
tially) complete suppression of the metallic temperature 
dependence occurs at fields below the full spin polariza- 
tion field (at finite temperature) where the 2D system 
is still partially spin-polarized. In Fig. 3 wc show the 
calculated spin-polarization as a function of the in-plane 
magnetic field for different temperatures. As the tem- 
perature increases the net spin polarization at B = Bg 
decreases due to the thermal excitation. The suppres- 
sion of "metallicity" in this simple explanation arises 
from two competing temperature- induced mechanisms in 
the partially spin-polarized 2D system at parallel fields 
just below complete saturation. The two mechanisms 
counter-balance each other because one increases screen- 
ing and the other decreases screening. Below (but close 
to) the full spin-polarization field B < Bg, increasing 
temperature reduces screening through the direct thermal 
broadening, but enhances screening by thermally exciting 
reversed-spin quasiparticles (i.e. by thermally reducing 
the net spin polarization). Actually, these two compet- 
ing temperature-induced mechanisms are always present 
in the 2D carrier system at any finite parallel field: in- 
creasing temperature induces two competing trends in 
screening - increased thermal broadening reduces screen- 
ing and thermal excitation enhances screening by reduc- 
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as the spin degeneracy. In (a) the lines correspond to the re- 
sults for T = 0.0, 0.5. 1.0 2. OK (from bottom to top). In (b) 
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show the normalized resistivity as a function of temperature 
at fixed magnetic fields B« =0, Bs, and 2Bs. 



ing the net spin polarization. Our theory and the numer- 
ical results, of course, include these two mechanisms. 

In Fig. 4 we show the calculated polarizabilities, 
IL(q,T)/NFO (Npo = 2Nf = .grm/27r), for different in- 
plane magnetic fields; (a) for _B = and B — 2Bs, and 
(b) for B = Bg. At B = (i.e., the system is unpolar- 
ized and n+ = n_), the decrease of the screening func- 
tion (i.e. n) with increasing temperature at q = 2kpo 
{kpo = \j2TmJgy) gives rise to increasing resistivity as 
the temperature increases. Most of the temperature de- 
pendence of the resistivity at low temperature regime 
comes from the suppression of screening &i q = 2kpo. 
The reduction of screening near q — 2kpo determines the 
strength of metallicity at i? = 0. At B = 2Bg (i.e., the 
system is completely polarized and n+ = n at T = 0), 
we find a very similar behavior of the screening func- 
tion, but the rate of the change of screening at q = 2kp+ 
{kp+ = ^JA^-Knjg^ is reduced by about a factor of four 
compared to the B = case in the same temperature 
range. Thus, we expect that the temperature depen- 
dence of /9(T) at i? = is roughly a factor of 4 stronger 
than that at S = 2Bs. At B = 2Bs the population of the 
minority band due to the thermal excitation is negligible 
in the given temperature range. In Fig. 4(b) we show the 
polarizability function at, B = Bg. As the temperature 
increases the screening function decreases at q = 2kp+ 



for B = Bs, but due to the occupation of the minority 
band we also find an enhancement of the screening func- 
tion near q = 2kp- {kp- = y/Awn~fg^ < kpo). Thus, 
the enhancement of the screening function near 2qkp- 
with increasing of the minority spin carrier density and 
the reduction of screening near q = 2kp+ by thermal 
broadening give rise to weak temperature dependence in 
p around B ^ Bg. 

It is, in principle, possible for these two competing 
mechanisms to almost completely cancel each other, close 
to B Bg. In Fig. 5 we show our numerical results for 
both situations of the valley degeneracy being affected 
(and unaffected) by the parallel field. Obviously the tem- 
perature independence of p{T) aromid B ^ Bg is more 
prominent in the situation where the valley degeneracy 
(as well as the spin degeneracy) is lifted by the external 
magnetic field. We emphasize that these competing tem- 
perature induced trends in screening in the presence of a 
finite spin-polarizing parallel magnetic field (i.e. the di- 
rect suppression of screening by increasing temperature 
and the indirect enhancement of screening by the tem- 
perature induced decrease of effective spin-polarization) 
is a generic phenomenon, and is therefore present in all 
2D carrier systems. In particular, these competing mech- 
anisms could be the underlying cause for the occurrence 
of the so-called parallel field-induced 2D metal-insulator 
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transition, which has been widely discussed in the context 
of hole transport in high-mobility 2D p-GaAs structures. 
It may, therefore, be possible that the peculiar temper- 
ature independence of 2D resistivity around B ^ Bs 
observed in n-Si MOSFETs arises from both the par- 
allel field induced lifting of valley degeneracy as well as 
the competing mechanisms discussed above (since valley 
degeneracy is the only physical mechanism that distin- 
guishes the Si MOS system from all the other 2D semi- 
conductor systems). 

It is important to emphasize that our finding (see Fig. 
5) of the temperature dependence of p(T, B) being the 
weakest below the full-spin-polarization field (at finite 
temperature), i.e. for B m Bs-, and the strongest for 
B = Q with the fully spin-polarized system {B Bs) 
having an intermediate behavior, as observed experimen- 
tally in ref. [23], is valid independent of whether the Si 
valley degeneracy is lifted by the applied magnetic field 
or not. This is a generic result which is always true in all 
2D systems independent of the number of valleys; The 
temperature dependence of p{T, B) is always the weakest 
for B K Bg due to the competition between the thermal 
excitation between spin up/down bands and the ther- 
mal suppression of screening. This is obvious from Fig. 
5 where the temperature dependence of p{T^ B) is the 
weakest around B ^ Bs, both in Fig. 5(a) and in Fig. 
5(b) with the assumed valley degeneracy being 2 and 
1 respectively. As we have emphasized above, the ad- 
ditional assumption of a magnetic field induced lifting 
of the valley degeneracy makes p{T) almost a constant, 
completely independent of temperature, in quantitative 
agreement with the data of Tsui et al. [7] . Whether this 
actually happens or not can only be decided experimen- 
tally [22]. We should, however, emphasize that even if 
the valley splitting A„ is small (e.g. A„ < 2Ef so that 
the valley degeneracy is not completely lifted), a partial 
lifting of valley degeneracy (A^, ^ 0) will certainly con- 
tribute to the physics being proposed here. For a partial 
lifting of the valley degeneracy, our calculated temper- 
ature dependence of p{T) will be intermediate between 
that shown in Figs. 1(b) and 1(c) or Figs. 5(a) and 5(b) 
as the case may be. In particular, a finite valley split- 
ting will lead to the same competing thermal trends in 
p{T) as discussed above and by Shashkin et al. [23] for 
the spin splitting, producing an additional suppression of 
p{T) bringing our theoretical results to better quantita- 
tive agreement to the Tsui at al. [7] data. 

IV. DISCUSSION 

We have shown in this paper that the enigmatic tem- 
perature dependence of the parallel field magnetoresis- 
tance, p{T,B), in 2D Si MOSFETs can be qualitatively 
and semiquantitatively well explained by the screening 
theory, assuming charged impurity scattering to be the 



dominant low-temperature resistive scattering mecha- 
nism. Our screening theory based explanation, in fact, 
becomes quantitatively accurate if we assume that the 
valley degeneracy is lifted by the external magnetic field 
in the same manner as the spin degeneracy. The fully 
spin- and valley-polarized Si MOS system at high field 
has a factor of 4 lower density of states than the corre- 
sponding zero-field spin- and valley-unpolarized system, 
leading consequently to substantially weaker screening 
with the corresponding strong suppression in the temper- 
ature dependence of the resistivity. We show that, inde- 
pendent of the valley degeneracy question, the weakest 
temperature dependence in the resistivity happens not 
dX B ^ Bg (where the carriers are essentially completely 
spin-polarized), but at B < Bg, just below the T = 
full spin-polarization field, where the system is almost, 
but not quite, fully spin-polarized. Thus, the metallicity 
(i.e. the magnitude of dp/dT) is the strongest for B = 
and the weakest for B < Bg with the fully spin-polarized 
{B ^ Bg) situation being intermediate. This specific 
prediction of the screening theory is in good agreement 
with experimental observations both in n-Si MOSFETs 
[23] and in n-Si/SiGe [20] 2D systems. The key feature of 
our theory is a realistic and quantitatively accurate de- 
scription of 2D screening taking into account finite tem- 
perature and finite spin-polarization effects on an equal 
(and non-perturbative) footing. The important ingredi- 
ent of physics leading to the strong suppression of metal- 
lic temperature dependence of p{T,B « Bg) in our the- 
ory is a competition between two physical mechanisms: 
Thermal excitation of reversed spin quasiparticles lead- 
ing to enhanced screening with increasing temperature 
for B < Bg and the direct thermal suppression of screen- 
ing by the majority spin carriers. We note that this com- 
petition would always lead to a strong suppression in the 
temperature dependence of p{T) around B < Bg, inde- 
pendent of the member of valleys involved in the 2D sys- 
tem. It is also important to emphasize that for B > Bg 
(and at not too low carrier densities) our theory pre- 
dicts a negative magnctorcsistancc with monotonically 
decreasing p{B) with increasing B due to the Coulomb 
matrix element effect arising from an increasing effective 
kF{B) in the system. This also leads to the suppression 
of metallic temperature dependence at high fields. 

Our finding of the qualitative agreement between 
the screening theory and the experimentally observed 
p{T\ B) in Si MOSFETs is of considerable importance. 
Much has been made in the literature of the claimed non- 
metallic (i.e. temperature independence or even insulat- 
ing) behavior of 2D magnetoresistatce in Si MOSFETs at 
high {B > Bg) parallel fields. In particular, interaction- 
based theories [15] of 2D metallicity have been claimed 
[1] to provide the definitive explanation for 2D metallicity 
in Si MOSFETs because the interaction theory predicts 
the fully spin and valley-polarized 2D system to have 
an insulating (with p oc l/T) temperature dependence, 
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and it has been categorically asserted that the screen- 
ing theory cannot possibly be qualitatively correct [24] 
since it predicts a (weak) metallic behavior in the fully 
spin-polarized ballistic regime. We have shown in this 
work that the screening theory is, in fact, in good quali- 
tative (and even quantitative, if we allow the possibility 
of valley degeneracy lifting) agreement with the most re- 
cent Si MOS experimental data [7,23] in the presence of a 
parallel magnetic field, and in particular, consistent with 
our theoretical predictions, the strongest suppression of 
the metallic temperature dependence of resistivity hap- 
pens for B < Bg where the 2D system is almost, but 
not quite, fully spin polarized. Thus, p{T, B) manifests 
the strongest metallicity at S = (when the system is 
unpolarized) the weakest aX B Bg (when the system 
is almost polarized), and intermediate at B ^ Bs (when 
the system is fully polarized) both in the screening the- 
ory and in the experiment. Exactly the same behavior 
has also been seen recently by two experimental groups 
[20] in 2D n-Si/SiGe electron systems which we have 
theoretically explained [25] recently using the screening 
theory. Thus, the parallel field induced suppression of 
metallic temperature dependence in the ballistic regime 
is now theoretically understood (at least qualitatively) 
for both n-Si MOS and n-Si/SiGe 2D electron systems, 
and the early experimental discrepancies between these 
two systems have now been resolved [23] with both sys- 
tems agreeing reasonably with the screening theory pre- 
dictions in the ballistic regime. 

There are various proposed non-Fermi liquid scenarios 
for the 2D metallicity which depend crucially on impor- 
tant distinctions between the spin polarized {B > Bg) 
and the unpolarized {B = 0) transport regimes. At least 
one of these non- Fermi liquid scenarios, based on a spec- 
ulative coexistence between 2D Wigncr crystal and elec- 
tron liquid phases [26] is apparently invalidated by the 
observed weak metallicity in the spin-polarized Si MOS- 
FETs, since the theory [26] specifically predicts abso- 
lutely no temperature dependence in /Ci(T; B ^ Bg) \j\ 
the spin-polarized phase in direct contradiction with the 
recent experimental results. [20,23] 

We briefly discuss the role of interaction in our screen- 
ing theory, touching upon the closely related (and im- 
portant) question of why the Boltzmann-RPA screening 
theory seems to provide a good description of the 2D 
transport properties at low carrier densities where the 
dimensionless interaction parameter (oc nT^I'^ in 2D), 
defined as the ratio of the average Coulomb energy to the 
Fermi energy at T = 0, is larger than one (in Si-based 2D 
systems of interest in this paper, Tg ~ 6 in the experimen- 
tally relevant carrier density range of n ~ IQ-^-'^cm"^). 
We believe that one possible reason for the success of 
the RPA screening theory is that the dominant contri- 
bution to the 2D resistivity arises from charged impu- 
rity scattering, and RPA-Boltzmann theory does an ex- 
cellent job of regularizing the associated Coulomb dis- 



order through the self-consistent screening model. The 
2D transport problem is, in fact, semiclassical since the 
temperature, expressed in units of the Fermi tempera- 
ture, is not necessarily small (often T/Tp ^ 0(1)), again 
making finite temperature RPA an excellent approxima- 
tion. It is important to emphasize that at a fixed finite 
temperature T, decreasing 2D density increases both 
and T/Tp (since Tp cx: n), and the finite-temperature 
interaction parameter, e.g. rs/{T/Tp) ^ ^Jn, actually 
decreases with decreasing density! Therefore, it is not 
obvious at all that finite temperature 2D transport be- 
havior increasingly becomes strong-coupling like as elec- 
tron density is decreased at a fixed temperature since the 
low-density (r^ oo) limit is also at the same time the 
classical infinite temperature {T/Tp — > oo) limit. Under 
these circumstances, the semiclassical RPA-Boltzmann 
finite temperature theory may, in fact, become increas- 
ingly more valid as the carrier density is reduced, partic- 
ularly since the carrier temperature can often not be re- 
duced below 100 mK due to electron heating effects. This 
finite-temperature aspect of the low-density 2D transport 
problem, which has not been adequately emphasized in 
the literature, may very well be playing an important role 
in making RPA a particularly good approximation in the 
2D MIT phenomena. 

We note that RPA-based many-body theories also 
seem to describe reasonably well [27,28] the experimental 
behaviors of finite-temperature, dilute 2D carrier systems 
with respect to collective mode dispersion [29] and inter- 
layer drag measurements [30]. This finite-temperature 
(actually high-temperature) aspect of this problem most 
likely also invalidates any possible relevance of Wigner 
crystallization to the physics of 2D MIT phenomena - 
a Wigner crystal would thermally melt [27,31] at the 
"relatively" high temperatures (i.e. T/Tp ^ 0(1)) at 
which the low density 2D transport experiments are typ- 
ically carried out. It needs to be emphasized that the 
asymptotic low-tcmpcrature behavior of p{T) is essen- 
tially never observed in 2D transport experiments since 
the experimental p{T) always saturates at the lowest 
measurement temperatures (often for T < lOOmK) in- 
dicating that the electrons may not be cooling down to 
the T/Tp <^ 1 regime. In this intermediate to high 
temperature regime our screening theory, which is non- 
perturbative in temperature, may very well be the qual- 
itatively correct approximate theory. Understanding the 
true asymptotic temperature dependence of p{T) in the 
T/Tp <C 1 regime, where interaction effects must be im- 
portant in a low-density strong-coupling electron system, 
remains an important open experimental and theoretical 
challenge. 

We can speculate on the possibility of incorporating 
interaction effects into our screening theory of 2D trans- 
port. The 2D resistivity, p{n,T,B), depends on density, 
temperature, and the applied magnetic field, which may 
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be expressed in the dimensionless units to write 

p{n,T,B) = p{qTF/'2kF,T/TF,A,/TF), (11) 

where As = gpsB is the spin spHtting induced by the 
external magnetic field. We note that ^tfI'^^f oc m; 
T/Tf oc m; and Ag/Tp oc gm oc x where m, g, and x 
are respectively the 2D effective mass, Landc g-factor, 
and spin susceptibility. In the RPA theory, interaction 
effects are ignored in the effective mass and the ^r- factor, 
and therefore, the nonintcracting effective mass and the 
susceptibility are used in calculating the resistivity. In 
the spirit of the Landau Fermi liquid theory, we could 
crudely incorporate interaction effects in the theory by 
using the quasiparticle effective mass (m* ) and the quasi- 
particle spin susceptibility (x*) or equivalently the quasi- 
particle g-factor (g*) in the effective RPA calculation. A 
rigorous justification for such a Fermi liquid renormal- 
ization of the effective mass and g'-factor in the RPA- 
Boltzmann transport theory is unavailable, and indeed 
all our numerical results utilize the non-interacting band 
mass and g-factor in the Boltzmann-RPA calculations, 
but we speculate that such a physically motivated ap- 
proximation (i.e. TO* and g* replacing to and g) to incor- 
porate interaction effects may be quite reasonable since 
the observed or measured quantities are actually m* and 
g* (and x* = fn*g*) and not the bare quantities (m 
and g). It is interesting to point out that such an ad 
hoc approximation scheme incorporating interaction ef- 
fects (by using to* and g* in the Boltzmann-RPA trans- 
port theory) does considerably improve the quantitative 
agreement between theory and experiment. In particu- 
lar, the interacting susceptibility x* is renormalized [32] 
by a factor of 3 (i.e. x* Ix 3) at n ^ IO^^cto^^ in 
Si MOSFETs leading to the theoretical spin-polarization 
saturation magnetic field Bg (~ 3 — 4T) being in good 
agreement with the experimental data. Similarly, the use 
of the quasiparticle effective mass to* ^ 3to, at low densi- 
ties considerably increases the effective values of q^p/2kF 
and T/Tp (where "starred" quantities use to* rather than 
to) , again bringing theory and experiment in good quan- 
titative agreement. Whether such an ad hoc "improve- 
ment" of the theory using quasiparticle rather than bare 
Fermi liquid parameters can be rigorously theoretically 
justified or not remains an important open question for 
the future. We note, however, that the screening the- 
ory obtains semi-quantitative and qualitative agreement 
with the existing 2D transport data even without any 
such "renormalization" . 

Finally, we note that there still seems to be some quan- 
titative difference in the observed experimental temper- 
ature dependence of p{T; B) between Si MOS [7,8] and 
Si/SiGe [20] 2D electron systems. For example, in the Si 
MOS system p{B) is enhanced approximately by a fac- 
tor of three at low temperatures as B increases from zero 
to Bs whereas in the Si/SiGe system [20] the enhace- 
ment factor is only almost 1.8. This difference arises (at 



least partially) from the difference in the charged impu- 
rity distribution in the two systems — in the Si/SiGe 2D 
electron system the dominant scattering is from remote 
dopants and background impurities [25] whereas in Si 
MOS system the scattering is mostly by interface charged 
impurities and interface roughness. In addition, the finite 
magnetic field behavior in these two Si-based 2D systems 
may also differ by virtue of the valley degeneracy being 
lifted in the Si MOS system, but not in the Si/SiGe sys- 
tem. Such a non-universal lifting of the valley degeneracy 
is certainly possible since it is well-known that the Si val- 
ley degeneracy near an interface depends critically on the 
microscopic details of the interface, and it is entirely pos- 
sible for the valley degeneracy to be lifted in the Si MOS 
system, but not to be lifted in the Si/SiGe system since 
the latter system has a much better atomically smooth 
epitaxial interface. In fact, the non-universal lifting of 
valley degeneracy may apply even to different Si MOS 
2D systems, and could provide one possible underlying 
reason for the observed difference [7,8] in the tempera- 
ture dependence of p{T;B > Bg) in Si MOS data from 
different groups. Whether such a non-universal parallel 
field-induced electron valley degeneracy lifting actually 
happens in reality can only be decided experimentally; we 
are only suggesting here the theoretical possibility based 
on our analysis of the experimental transport properties 
oip{T;B). 

V. CONCLUSION 

In conclusion, we consider theoretically the parallel 
magnetic field induced suppression of the screening of 
long range bare Coulomb disorder in Si MOSFETs, show- 
ing that the experimentally observed strong suppression 
of 2D metallicity in the temperature dependence of the 
magnetoresistance can be qualitatively and semiquanti- 
tatively understood as arising from spin (and perhaps 
even valley) polarization induced reduction in carrier 
screening, leading to stronger parallel field-depc>nd(uit (ef- 
fective disorder in the system. The competing meclia- 
nism of direct thermal reduction of screening and indi- 
rect enhancement of screening through the thermal sup- 
pression of spin polarization may also be playing an im- 
portant quantitative role in the temperature dependence 
of p{T;B ^ Bs). Our theory, as presented in the cur- 
rent work and in our recent work [25] on the 2D Si/SiGe 
electron system along with the recent experimental work 
[20,23], may resolve the earlier discussed qualitative dis- 
agreement between Si MOS and Si/SiGe 2D systems, es- 
tablishing the same qualitative behavior in all 2D Si sys- 
tems. More work is still needed to precisely understand 
interaction effects on 2D transport by going beyond the 
physically motivated screening theory of our work. 
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